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year.  It  was  found  that  while  climatological  maps  oftfjT3' 
calculated  from  filtered  data  show  Spatially  coherent 
patterns  that  are  clearly  associated  VHUi^jM-nth-ly  storm 
tracks,  similar  maps  of  c'ufl'3'’do  not  show  such  coherent 
patterns.  V _ ^ 

Maps  of  anomalous  storm  activity,  as  given  by  the  dif¬ 
ference  between  the  monthly  value  and  the  climatological 
value  for  that  month,  are  also  prepared  and  qualitatively 
compared  to  observations  of  large  scale  sea  surface  tem¬ 
perature  (SST)  anomalies^  Fairly  good  agreement  is  found 
between  monthly  mean  u*3  i^nomal ies  calculated  from  the 
filtered  data,  and  the  observed  development  of  SST  anoma¬ 
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ABSTRACT 


Synoptic  storm  activity  over  the  North  Pacific  Ocean  is 
investigated  using  a  special  series  of  six-hourly  surface 
wind  analyses  prepared  by  Fleet  Numerical  Oceanographic  Cen¬ 
ter  (FNOC)  covering  the  period  1969-1973.  Temporal  variance 
of  the  surface  wind  components  at  each  grid  point  is  "high 
pass"  filtered,  thereby  removing  all  temporal  variability 
except  that  having  time  scales  less  than  ten  days.  Both 
the  total  and  the  filtered  wind  components  were  used  to 
calculate  the  cube  of  the  friction  velocity  (u*  )  and  the 
wind  stress  curl  (curl  t). 

Monthly  "climatologies"  are  calculated  by  computing  ten- 
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year  averages  of  monthly  mean  u*  and  curlzt  from  total  and 
filtered  wind  components  for  each  month  of  the  year.  It  was 

3 

found  that  while  climatological  maps  of  u*  calculated  from 
filtered  data  show  spatially  coherent  patterns  that  are 
clearly  associated  with  monthly  storm  tracks,  similar  maps 
of  curlzx  do  not  show  such  coherent  patterns. 

Maps  of  anomalous  storm  activity,  as  given  by  the  dif¬ 
ference  between  the  monthly  value  and  the  climatological 
value  for  that  month,  are  also  prepared  and  qualitatively 
compared  to  observations  of  large  scale  sea  surface  tem¬ 
perature  (SST)  anomalies.  Fairly  good  agreement  is  found 

3 

between  monthly  mean  u*  anomalies  calculated  from  the 
filtered  data,  and  the  observed  development  of  SST  anoma¬ 
lies  during  the  winter  of  1976-1977. 
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I.  INTRODUCTION 


There  have  been  many  studies  investigating  the  response 
of  the  upper  ocean  to  strong  atmospheric  forcing  events.  By 
using  many  years  of  data  from  Ocean  Weather  Ships  P,  V,  and 
N,  Elsberry  and  Camp  (1973)  showed  that  a  major  fraction  of 
the  oceanic  thermal  structure  modification  during  September 
to  December  is  due  to  strong  atmospheric  forcing.  An  earlier 
study  by  Simpson  (1969)  suggests  that  synoptic  scale  distur¬ 
bances  control  air-sea  interaction  over  the  extratropical 
oceans.  The  choice  of  an  indicator  of  the  atmospheric  circu¬ 
lation  has  varied.  Namias  (1972  and  1978)  worked  with 
monthly  mean  700  mb  height  fields  while  Davis  (1976)  worked 
with  monthly  mean  sea-level  pressure  fields.  In  the  latter 
study  Davis  used  linear  statistical  predictors  to  examine 
connections  between  thermal  variability  in  the  North  Pacific 
Ocean  and  variations  in  the  state  of  the  overlying  atmosphere. 
He  also  suggests  that  atmospheric  storminess  may  show  a 
greater  connection  with  sea  surface  temperature  (SST)  than 
does  pressure.  In  a  preceding  thesis  study,  using  12-hourly 
data  covering  but  a  single  year,  Heise  (1977)  found  that 
synoptic  storms,  or  the  lack  of  them,  were  related  to 
anomalously  cold,  or  warm,  SST  patterns.  Heise  was  not, 
however,  able  to  examine  anomal ous  storm  activity  nor  to 
compare  it  with  SST  anomalies. 
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The  objective  of  this  study  has  been  to  examine  ten 
years  of  synoptic  wind  data  for  the  purpose  of  describing 
synoptic  storm  activity.  Using  a  series  of  six-hourly  sur¬ 
face  wind  analyses  prepared  by  Fleet  Numerical  Oceanographic 
Center  (FNOC)  covering  the  period  1969-1978,  the  temporal 
variance  of  the  surface  wind  components  at  each  individual 
grid  point  was  "high-pass"  filtered.  The  filter  removed 
most  of  the  temporal  variability  having  time  scales  greater 
than  ten  days.  Therefore,  the  circulation  primarily  due  to 
moving  synoptic  storm  systems  was  extracted  from  the  total 
circulation. 

Heise  (1977)  showed  that  filtered  data  does  describe 
monthly  mean  storminess  over  a  period  of  one  year  (1975)  and 
that  this  storminess  was  qualitatively  related  to  SST 
anomalies  during  that  year.  From  both  the  filtered  and  un¬ 
filtered  wind  data,  Heise  (1977)  computed  the  vertical  compo¬ 
nent  of  the  curl  of  the  wind  stress  (curlzx)  and  the  friction 

3 

velocity  cubed  (u*  ).  These  quantities  were  selected  because 
of  the  well-established  relationship  between  curlzx  and  Eknan 
convergence  and  divergence,  and  the  relationship  between 

3 

u*  and  the  mechanical  mixing  of  the  sea  by  the  overlying 
atmosphere.  This  same  analysis  was  used  in  this  study  to 
identify  storm  tracks  and  to  determine  the  extent  to  which 
synoptic  storm  activity  mixes  and/or  pumps  the  ocean. 

From  the  ten  years  of  six-hourly  surface  wind  data,  a 

3 

monthly  climatology  of  both  u*  and  curlzx  was  established. 
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Using  this  cl  imatol  ogy,  u*  and  curl^  anomalies  were  also 
calculated  for  each  month  of  the  ten  year  period.  The  re¬ 
mainder  of  this  paper  is  devoted  to  describing  the  proce¬ 
dures  used  in  data  processing,  data  analyses,  and  interpret¬ 
ing  the  patterns  of  climatological  u  *  ^  and  curlzx.  A  case 
study,  which  compares  the  change  in  anomalous  SST  during 

3 

the  winter  season  1976-1977  to  the  averaged  anomalous  u* 
and  cur1zx  for  that  period,  will  also  be  discussed. 


II.  DATA  SET  DESCRIPTION  AND  ANALYSES 


A.  DATA  ACQUISITION  AND  PROCESSING 

The  wind  data  were  obtained  from  FNOC  through  Scripps  In¬ 
stitution  of  Oceanography  and  were  based  on  six-hourly  surface 
pressure  analyses  over  the  North  Pacific  Ocean.  The  first 
seven  years  (1969-1975)  of  data  were  produced  by  a  private 
company  under  a  contract  from  FNOC  (Mendenhall,  Holl  and 
Cuming,  1978),  with  the  last  three  years  (1976-1978)  coming 
directly  from  FNOC  operational  analyses.  The  surface 
pressure  analyses  were  derived  by  using  previously  computed 
surface  pressure  analyses  as  a  first  guess  and  upgrading 
the  analyses  with  surface  pressure  and  wind  velocity  obser¬ 
vations  from  the  U.S.  National  Climatic  Center  land  and  ship 
data.  The  winds  were  adjusted  to  an  elevation  of  19.5 
meters  above  sea  level.  However,  one  weakness  in  using 
surface  pressure  analyses  is  that  an  underestimation  of 
high  wind  speeds  may  result  (Lazanoff  and  Stevenson,  1978). 
From  these  data,  u  (eastward)  and  v  (northward)  wind  com¬ 
ponents  were  extracted  from  grid  points  on  the  FNOC  63x63 
Northern  Hemisphere  polar  stereographic  array  between  0°  and 
70°N  and  from  120°W  to  100°E  (North  Pacific  Ocean  area). 

The  grid  spacing  varies  from  381  km  at  60°N  to  208  km  at  the 
equator.  At  each  grid  point  there  were  four  daily  observa¬ 
tions  (003MT,  06GMT,  1 2  GMT ,  18GMT)  throughout  the  ten-year 
peri od . 
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B.  TREATMENT  OF  MISSING  DATA 

The  maximum  number  of  consecutively  missing  observations 
was  seven  and  this  happened  twice  in  the  ten-year  period. 
Overall,  less  than  six  percent  of  the  total  observations  was 
missing.  Figure  1  shows  the  percent  of  missing  data  each 
month  for  the  ten  year-period.  It  should  be  quite  apparent 
from  the  figure  when  the  special  analysis  produced  by  the 
private  company  terminated  (.1  975  )  and  when  the  operational 
data  from  FNOC  began  (1976).  All  missing  data  groups  were 
filled  by  interpolating  linearly  from  preceding  and  succeed¬ 
ing  map  times. 

C.  DETRENDING  AND  FILTERING  OF  DATA 

The  next  step  was  to  detrend  and  filter  the  data.  By 
detrending  the  data  any  linear  trends  were  removed.  After 
the  appropriate  slopes  and  intercepts  were  calculated,  the 
equation  used  in  detrending  was 

Yn,  ’  *  b>  * 1 1-1  * 

where  Ym  is  the  detrended  value  at  the  mth  time  observation, 
ym  is  the  original  observation,  a  is  the  slope,  b  is  the 
intercept  and  t  =  (m-l)At,  with  At  =  six  hours,  being  the 
time  interval.  The  method  of  least  squares  applied  to  each 
grid  point  forms  the  basis  of  the  detrending  program  used. 

In  order  to  extract  synoptic  storm  activity  from  the 
ten-year  period,  a  41-point  symmetric  high-pass  filter  was 
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applied  to  the  detrended  time  series  of  velocity  components 
at  each  grid  point.  The  filtering  equation  was 


20 

ui,j^tm^  “  woui,j^tm^  +  ^ u  i  ,  j  ^  tm+k  ^ +u  i  ,  j  ^  tm-k  ^  wl 

k  =  l 


( 1 1-2 ) 
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k  =  l 


where  u.  .(t  )  and  v.  .( tm)  represent  the  u  and  v  wind  com- 

ponents  respectively  at  the  grid  point  (i,j)  at  the  time 

t  ,  and  w,  are  the  coefficients  for  the  filter.  The  coeffi 
m  k 

cients,  shown  in  Table  I,  were  computed  using  the  following 
Gaussian  filter, 


wQ  =  1  -At(  2ttc2  ) 

wk  =  -At(2iT  o2)"1*  exp(-tk2/2o2)  ,  ( 1 1  -  3 ) 

where  tk  =  kAt,  k=l,...20,  a  =  1.5  days  and  At  =  .25  days. 
The  response  of  the  filter  is  shown  in  Figure  2.  The  filter 
is  sensitive  to  periods  of  ten  days  or  less  or  frequencies 
greater  than  0.10  cycles  day’\  The  significant  response 
is  for  periods  of  less  than  five  days.  This  is  designed  to 
extract  events  within  the  synoptic  time  scale. 


TABLE  I.  Filter  coefficients  and  corresponding  weights, 
k  wk 


0 

0.933510 

1 

-0.065573 

2 

-0.062897 

3 

-0.058678 

4 

-0.053241 

5 

-0.046985 

6 

-0.040328 

7 

-0.033666 

8 

-0.027335 

9 

-0.021586 

10 

-0.016580 

11 

-0. 01 2385 

12 

-0.008998 

13 

-0.006359 

14 

-0.004370 

15 

-0.002921 

16 

-0.001899 

17 

-0.001201 

18 

-0.000739 

19 

-0.000442 

20 

-0.000257 

III.  CALCULATIONS  PERFORMED 


A.  EQUATIONS 

3 

Calculations  of  friction  velocity  cubed  (u*  )  and  the 
curl  of  the  wind  stress  (curlzx)  were  made  four  times  a 
day  at  every  grid  point  using  both  the  total  and  filtered 
wind  components.  The  equations  employed  for  the  calcula¬ 
tions  were: 

u*3  =  £ C D ( u 2  +  v2)  ]3/2  ( 1 1 1 -1  ) 

curV  ■  "2(-k  ^  -  ly  's25-5  ]  (III-2) 


where , 


(  0.75  +  0.00067(u2  +  v2)J'2)  10‘3 
( 1  +s  i  n<p  t )  /  (1  +s  i  n<p ) 

-x,  sinA*  -  x  cosA* 

A  *  ( II 1-3) 

x,  cosA*  -  x  sinA* 

PaCD(u2  +  v2)i5u 

PaCD(u2  +  v2)!sv 

XT / 2  -  2  arctan[(x2+y2)li/a(l+sin$t)  ] 
arctan(y/x)  . 


In  the  above,  a  is  the  radius  of  the  earth,  $>t  is  the  true 
latitude  of  the  polar  stereographic  projection,  60°,  and 
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■in  i 


-3 

the  air  density  p  was  taken  as  a  constant  1.3  x  10  gm  cm 

a 

Other  expressions  used  included  the  non-dimensional  drag 
coefficient  CQ  (Garratt,  1977),  the  zonal  and  meridional 
stress  components  x.  and  t  ,  and  the  map  factor  m,  which 

A  $ 

represents  the  coordinate  tra ns  formation  from  spherical  co¬ 
ordinates  (A,?)  to  polar  stereographic  coordinates  (x,y) 

(see  Haltiner,  1971,  p.13).  The  same  equations  were  used 
on  both  the  filtered  and  total  wind  fields  with  the  appro¬ 
priate  u  and  v  components.  In  the  finite  difference  approxi 
mations  used  to  calculate  curlzx,  centered  differences  on  a 
staggered  grid  on  the  polar  stereographic  projection  were 

used.  If  xY  and  xv  are  defined  on  the  array  i=l,... 
i.j  i,J 

IM  and  j=l,...JM,  then  curlzT  is  defined  on  a  staggered  grid 
by 


2,1, 


(curV>1,J  =  l2<(,Y/m)Ul.j*l  + 

-  jdtx/in),  J.i  +  ,j*l) 


where  m,.  j  is  an  averaged  map  factor  defined  by 


»  .  ISiiLll!  -tfl.J*!)* 

mi,j  4 

and  where  the  range  of  i.j  is 


i =1 , . . . IM- 1 


( 1 1 1  -  5 ) 

and  j-1 , . . . JM-1 . 
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B.  u  3  AND  CURL  t  DEVELOPMENT 
z 

As  will  be  shown  later,  a  measure  of  the  monthly  mean 
storminess  (which  will  be  used  to  develop  a  ten  year  clima¬ 
tology)  and  the  corresponding  storm  tracks  can  be  calcu¬ 
lated  from  the  high-pass  filtered  wind  components.  Since 

3 

u*  represents  friction  velocity,  u*  was  used  because  it 

is  a  measure  of  the  mechanical  mixing  in  the  ocean  due  to 

the  overlying  atmospheric  circulation.  Elsberry  and  Camp 

3 

(1973)  obtained  a  high  correlation  between  u*  and  mixing 

3 

of  the  ocean.  A  u*  value  was  calculated  for  each  data 

point  and  for  each  observation.  For  comparison  purposes, 

3 

u*  calculations  were  performed  using  both  the  total  and 
the  filtered  wind  components. 

Curl  calculations  were  also  performed  on  both  the  total 
and  the  filtered  wind  fields.  The  objective  of  calculating 
curlzi  from  the  filtered  wind  data  was  to  see  if  there  was 
a  contribution  to  Ekman  pumping  related  solely  to  the 
monthly  mean  location  of  synoptic  storm  activity.  It 

3 

should  be  noted  that  the  above  calculations  of  u*  and 
curlzx  using  the  high-pass  filtered  wind  components  reveal 
the  effects  of  the  synoptic  activity  alone,  i.e.,  in  the 
absence  of  any  low  frequency  or  mean  flow.  None  of  the 
interaction  terms  (between  high  and  low  frequencies)  which 

3 

appear  in  the  nonlinear  quantities  u*  and  curlzx  are 
included.  This  point  is  discussed  in  the  conclusion  section 
below. 
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C.  DESCRIPTION  OF  CLIMATOLOGY  AND  ANOMALIES 

3 

With  u*  and  curlzx  values  calculated  for  each  grid 
point  and  for  each  observation,  a  monthly  mean  for  each 
month  between  January  1969  and  December  1978  was  calculated 
This  was  done  using  both  the  filtered  and  total  wind  compo- 

3 

nents.  It  should  be  noted  that  the  monthly  mean  u*  and 
curlzx  values  for  January  1969  and  December  1978  calculated 
from  the  filtered  data  consist  of  only  26  days.  This  was 
necessary  due  to  the  41-point  symmetric  high-pass  filter 
which  needed  5  days  (20  observations  in  time)  on  either 
side  of  the  point  to  be  filtered. 

With  the  monthly  means  computed,  the  next  step  was  to 
develop  a  climatology  and  then  to  use  the  climatology  and 
monthly  means  in  establishing  anomalies.  A  climatology 
based  on  the  ten  years  of  data  was  developed  for  each  month 

3 

of  the  year  by  averaging  all  the  u*  and  curlzT  monthly 
means  for  January,  then  February,  and  continuing  through 
December.  With  monthly  means  and  climatology  values,  anoma 

3 

lies  of  u*  and  curlzx  were  created  by  subtracting  the 
appropriate  climatology  from  each  monthly  mean  throughout 
the  ten-year  period.  In  the  next  section  attention  will  be 

3 

focused  on  climatological  values  of  u*  and  curlzx  as 
calculated  from  total  and  filtered  wind  components.  After 
that,  a  relationship  between  anomalous  storm  activity  and 
the  development  of  anomalous  SST's  will  be  examined. 
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IV.  DISCUSSION  AND  RESULTS 


A.  EXAMINATION  OF  u*3  and  CURL CLIMATOLOGY 

z 

3 

In  this  section,  values  of  climatological  u  *  calculated 
from  the  high-pass  filtered  wind  components  will  be  shown 
to  be  a  good  representation  of  storminess  and  storm  tracks 
over  the  ten-year  period.  On  the  other  hand,  values  of 
climatological  curlzT  are  not  as  conclusive.  The  u*  and 
curlzx  parameters  were  used  because  of  their  role  in  mixing 
and  pumping  of  the  ocean,  respecti vel y .  The  12  months  of 

3 

climatological  u*  patterns  calculated  from  the  total  wind 
field  will  be  examined  first,  followed  next  by  an  examina- 

3 

tion  of  the  climatological  u*  patterns  calculated  from  the 
filtered  wind  field.  The  same  procedure  will  be  followed 
using  c  url zt  resul ts  . 

3  3 

1.  Relationship  Between  u+  and  u,  from  Filtered  Data 
Figures  3-14  are  plots  of  climatological  u*3  calcu¬ 
lated  from  the  total  wind  field.  For  January  (Fig.  3)  the 

3 

maximum  values  of  u*  over  water  are  centered  at  approxi¬ 
mately  55°N  and  1 5 5 0 E  and  45°N  and  150°W.  In  February  (Fig. 
4)  there  is  little  change  with  the  two  maxima  remaining  in 
nearly  the  same  area  and  with  the  same  magnitudes.  By  March 
(Fig.  5)  the  u*  maxima  begin  to  decrease  in  magnitude. 
This  trend  continues  throughout  the  spring  and  into  the 
summer  (Figs.  6-11).  Tropical  disturbances  begin  appearing 
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in  May  (Fig.  7)  and  are  present  until  October  (Fig.  12). 

During  October  and  November  (Figs.  12,  13),  the  values  of 

3  •  3 

u*  begin  to  increase  with  areas  of  maximum  u*  located 

particularly  in  the  northwest  Pacific  and  a  secondary 

maximum  located  in  the  Gulf  of  Alaska.  The  maximum  values 

3 

of  u*  occur  from  late  fall  through  early  spring  with 

3 

little  movement  north  or  south.  The  u*  values  resulting 
from  the  trade  winds  south  of  30°N  are  too  small  to  appear 
in  the  figures  because  of  the  contour  levels  selected. 

The  values  of  u*  over  the  land  masses  (Siberia  and  Alaska) 
are  of  little  interest  in  this  study  which  is  focused  on 
synoptic  storm  activity. 

3 

Figures  15-26  are  plots  of  climatological  u*  cal¬ 
culated  from  the  filtered  wind  components.  Notice  that 
the  contour  levels  are  generally  an  order  of  magnitude 
smaller  than  in  the  corresponding  Figs.  3-14.  The  filtered 
wind  components  represent  storminess  and  monthly  mean  values 

3 

of  u*  from  filtered  data  are  a  measure  of  the  monthly  mean 
mixing  of  the  ocean  by  synoptic  storms.  In  January  (Fig. 

15)  a  well  defined  symmetric  and  coherent  u*  maximum  lies 
at  45°N  and  168°E.  If  maximum  winds  associated  with  a  given 
storm  occur  near  the  center  of  the  storm,  then  the  axis  of 

3 

maximum  u*.  represents  the  mean  position  of  the  storms  (storm 
tracks)  for  that  month.  In  this  case,  an  associated  storm 
track  runs  east-northeast  from  this  center  as  indicated  by 
the  dashed  lines.  The  months  of  February  through  April 
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(Figs.  16-18)  also  show  well  defined  and  coherent  patterns 

3 

of  u*  .  The  storm  tracks  continue  to  stand  out.  Also, 

3 

the  maximum  u*  values  remain  nearly  constant  during  this 
period.  It  is  not  until  May  (Fig.  19)  that  the  magnitude 
lessens  appreciably.  This  trend  of  decreasing  magnitude 
continues  through  the  summer  until  September  (Fig.  23)  is 
reached.  After  September  the  magnitudes  begin  increasing 
once  again.  This  can  be  expected  as  synoptic  storms  become 
more  developed  and  strengthen  during  early  fall.  In  the 
months  of  May  through  November  (Figs.  19-25)  tropical 
storm  activity  is  very  prevalent.  In  addition  to  the  de- 

3 

creasing  magnitude  of  u*  during  the  summer,  there  is  a 
slight  northward  and  eastward  movement  as  synoptic  storm 
activity  weakens  and  moves  north.  The  plots  clearly  show 
that  the  main  center  of  storminess,  throughout  the  year, 
lies  in  the  northwest  Pacific  Ocean  region. 

Now  the  question  may  be  raised  as  to  what  part  of 
the  total  u*3  field  is  due  to  the  synoptic  scale  part  of 
the  frequency  spectrum  as  determined  by  u*  calculated  from 
the  high-pass  filtered  wind  data.  Table  II  shows,  for  each 

3 

month,  the  percent  of  the  total  monthly  mean  u*  ,  that  is 
due  to  synoptic  storm  activity.  These  values  were  obtained 
by  dividing  the  u*  value  computed  from  the  filtered  wind 

3 

components  by  the  appropriate  u*  value  computed  from  the 
total  wind  components.  For  instance,  the  January  value  was 
calculated  by  dividing  the  maximum  oceanic  value  in  Fig.  15 


3 

TABLE  II.  Monthly  importance  of  filtered  u*  relative  to 
the  total  u*  field  in  percent. 


Month  Percent  Importance 


Jan 

33 

Feb 

34 

Mar 

34 

Apr 

37 

May 

20 

Jun 

20 

J  u  1 

20 

Aug 

30 

Sep 

40 

Oct 

30 

Nov 

32 

Dec 

50 

(13  x  104  (cm/s)3)  by  the  corresponding  maximum  in  Fig.  3 
(4  x  10  (cm/s)  )  resulting  in  0.33  or  33  percent.  It  was 
found  that  the  overall  average  was  32  percent.  During  the 
stormier  months  of  September  through  April  the  average  was 
36  percent  with  a  single  monthly  high  of  50  percent  in 
December.  Elsberry  and  Camp  (1978)  found  significantly 

3 

greater  u*  values  related  to  synoptic  storms.  The  discrep¬ 
ancy  is  most  likely  due  to  the  omission  of  interactions 
between  high  and  low  frequencies  in  the  u*  values  calculated 
from  filtered  data  in  this  study.  That  is,  the  calculations 
from  the  high  pass  filtered  data  used  in  this  study  show 
the  effects  of  storms  in  the  absence  of  a  mean  (or  low  fre¬ 
quency)  flow,  whereas  the  Elsberry  and  Camp  study  shows  the 
effects  of  the  storms  in  the  presence  of  a  mean  flow. 

The  major  result  of  the  above  calculations  is  two¬ 
fold  and  can  be  summarized  as  follows.  First,  the 


3 

climatological  values  of  monthly  mean  u*  calculated  from 

the  high  pass  filtered  wind  components  clearly  reveal  large, 

spatially  coherent  areas  of  synoptic  storm  activity.  Second, 

3 

the  magnitude  of  the  u*  which  comes  from  the  synoptic  part 

of  the  spectrum  alone  is  a  sizable  fraction,  approximately 

3 

one-third,  of  the  total  u*  . 

2.  Relationship  Between  Curl  x  and  Curl  x 
from  Filtered  Data 

Figures  27-38  are  the  monthly  climatological  maps 
of  curlz"r  calculated  from  the  total  wind  field.  Areas  of 
positive  curl  are  shaded.  A  strong  west  to  east  band  of 
positive  cur1zx  exists  in  all  months  between  40°-60°N.  These 
bands  are  especially  strong  during  the  winter  months  of 
December  (Fig.  38)  through  March  (Figs.  27-29)  with  a  north 
to  south  narrowing  of  the  band  during  the  summer  months  of 
May  through  September  (Figs.  31-35).  The  areas  of  maximum 
positive  curl^  are  widespread  but  appear  to  be  centered  in 
two  main  areas.  One  is  centered  at  55°N  and  160°E  with  the 
other  centered  at  55°N  and  140°W.  These  areas  become  less 
organized  in  the  spring  and  summer  months  of  March  through 
October  (Figs.  29-36).  This  broad  region  of  positive  wind 
stress  curl  coincides  with  the  location  of  the  semi -permanent 
Aleutian  Low  and  the  known  location  of  mid-latitude  cyclonic 
a  c  t  i  v  i  ty  . 

The  main  area  of  negative  curl^  lies  immediately 
south  of  the  positive  cu  rl z~  area  in  the  general  location  of 
the  semi  -  permanent  subtropical  high  between  20°-40°N.  This 
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west  to  east  band  has  no  seasonal  preference  in  either  area 
or  magnitude. 

The  monthly  maps  of  curl^t  calculated  from  the  high- 

pass  filtered  data  and  averaged  over  the  ten  years  are  shown 

in  Figs.  39-50.  A  broad  area  of  positive  cu  rl z?  is  found 

3 

in  mid-latitudes  similar  to  the  regions  of  filtered  u* 

shown  previously.  However,  the  patterns  are  much  more 

3 

"noisy"  and  somewhat  less  coherent  than  the  u*  patterns. 
During  the  months  of  January  and  February  (Figs.  39,  40)  the 
maximum  positive  curl^  areas  occur  at  40°N  between  150°- 
180°E.  As  spring,  represented  by  March  through  June, 

(Figs.  41-44)  approaches,  the  maximum  positive  curl^  shifts 
to  the  north  with  the  general  migration  of  the  mean  storm 
tracks.  The  curl^  reaches  its  maximum  positive  value  in 
March  (Fig.  41).  Beginning  in  April  and  continuing  through 
November  (Figs.  42-49)  the  center  of  maximum  cur1zT  remains 
quasi-stat ionary  between  40°-45°N  and  150°-170°E  before 
shifting  south  and  intensifying  in  December  (Fig.  50).  Also 
note  that  during  the  typhoon  and  tropical  storm  season  of 
June  through  October  (Figs.  44-48),  large  positive  and  nega¬ 
tive  areas  of  curlzT  exist  in  the  tropical  areas  especially 
off  Baja,  California  (10°-20°N,  100°-140°W)  and  in  the 
Philippine  and  North  China  Sea  areas  (10°-30°N,  120°-160°E). 

The  broad  band  of  negative  curlzx  to  the  north  of 
the  main  storm  track  behaves  in  the  same  manner  as  the  posi¬ 
tive  curl^.  During  the  months  of  December  (Fig.  50)  and 
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January  through  February  (Figs.  39,  40),  the  center  of 
the  region  of  negative  cur^t  lies  between  45°-55°N  and 
between  150°E-1F0°W.  In  the  summer  months  (Figs.  43-48)  it 
moves  north,  while  in  July  and  August  (Figs.  45,  46)  the 
negative  curl^  is  disorganized  and  very  weak.  Finally  the 
negative  curl^  becomes  more  organized  during  the  fall  (Figs. 
47-49)  . 

3 

As  with  u*  it  is  important  to  calculate  the  percent 
of  the  total  curlz^  field  due  to  storminess,  or  those  curlzT 
values  calculated  from  filtered  wind  components.  Table  III 
gives  a  monthly  percent  obtained  by  dividing  the  values  of 
curlzT  calculated  from  the  filtered  data  by  the  values  of 
curlzt  calculated  from  the  total  wind  field.  The  procedure 
used  in  obtaining  Table  III  was  the  same  as  used  for  Table  II. 

TABLE  III.  Monthly  importance  of  filtered  Curl zt  relative 
to  the  total  Curlzx  field  in  percent. 


Month 


Percent  Importance 


Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 


2 

3 

16 
1  3 
1  0 
1  5 
1  0 
10 
10 
1  0 
5 
3 


The  mean  value  for  the  entire  year  is  3.5  percent  which  is 

3 

significantly  lower  than  that  obtained  for  u*  . 


Further  comparisons  between  u*^  and  curlzT  calculated 

from  the  high-pass  filtered  data  show  a  major  difference  in 

3 

spatial  patterns.  The  u*  field  is  somewhat  more  coherent  and 

far  less  "noisy".  Synoptic  storm  areas  and  associated  tracks 

are  more  easily  discernible  and  the  patterns  are  very  smooth. 

3 

Since  storminess  contributes  relatively  more  to  the  total  u* 
field,  this  result  correlates  well.  The  major  result  of  the 
foregoing  examination  shows  that  the  contribution  of  synoptic 
storm  activity  to  vertical  turbulent  mixing  (u*  )  of  the 
ocean  is  significant,  while  the  contribution  of  synoptic 
activity  to  Ekman  pumping  (curlzT)  of  the  ocean  is  much  less 
significant.  In  the  next  section  an  example  of  SST  anomalies 

3 

will  be  compared  to  filtered  u*  and  filtered  curl^  anomalies. 

B.  CASE  STUDY  OF  ANOMALOUS  STORMINESS  AND  ANOMALOUS  SST 
It  is  reasonable  to  expect  that  at  the  same  location 
where  there  is  a  maximum  of  wind  mixing,  the  SST  decreases 
will  be  largest.  This  hypothesis  is  based  on  a  correlation 
of  upward  heat  flux  with  wind  speed  and  that  entrainment 
mixing  contributes  to  SST  decreases.  White  and  Haney  (1977) 
showed  that  during  the  autumn  and  winter  of  1976-77  the  storm 
forcing  and  SST  evolution  were  both  abnormal.  Between  30°-50°N 
in  the  central  portion  of  the  ocean,  the  SST  was  2  degrees 
colder  than  normal  during  the  period  in  which  the  frequency 
of  storms  was  4-6  times  greater  than  normal.  To  illustrate 
the  relationship  between  the  development  of  SST  anomalies  and 
anomalous  storm  activity,  the  SST  anomalies  that  developed 


3 

during  the  winter  season  of  1976-77  will  be  compared  to  u* 
and  curlzx  anomalies  calculated  from  the  filtered  wind 
components . 

1 .  Relationship  to  be  Tested 

3 

If  the  climatological  u*  or  curlzT  represents  the 
normal  or  average  synoptic  storm  activity  over  the  ten-year 
period,  then  departures  of  these  quantities  from  climatology 
may  be  related  to  SST  anomalies  through  anomalous  mixing 
and/or  pumping  of  the  ocean.  The  six-month  period  from  mid- 
September  1976  through  mid-March  1977  was  examined  to  test 
the  following  equation: 


(SST)  =  X  ,  (IV-1) 

where  SST  is  the  monthly  mean  sea-surface  temperature  anomaly 
(data  kindly  provided  by  Dr.  J.  Namias)  and  X  reDresents  the 

3 

monthly  mean  u*  or  curlzt  anomaly  calculated  from  the  fil¬ 
tered  data.  Integrating  (IV-1)  from  the  mid-September  1976 
(SEP76)  to  mid-March  1977  (MAR77),  the  equation  to  test  becomes 

^SST^MAR77  "  (sst)sEP76  =  *  (IV-2) 

where  J  is  the  integral  of  X  over  the  six-month  period.  This 
integral  was  calculated  from 

[FEB  77  -| 

l^XSEP76+XMAR77^  +  S  X  *  (IV-3) 

0CT76  J 

Using  (IV-2),  the  time  lag  between  atmospheric  forcing  and 
ocean  response  is  taken  into  consideration. 
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2.  Results  of  SST,  u  +  and  Curl.T  Anomalies 

Figure  51  (note  difference  in  latitudinal  coverage) 
is  a  plot  of  the  SST  anomaly  change  during  the  six-month 
period  (left  side  of  ( I V  -  2 ) ) .  Negative  values  are  shaded 
indicating  the  development  of  colder  than  normal  SST's,  while 
positive  values  are  unshaded.  The  main  area  of  cold  anomaly 
development  lies  between  25°-40°N  and  between  170°E-150°W. 

A  secondary  area  extends  southwestward  off  the  west  coast  of 
California  and  Baja.  Another  area  of  consequence  is  centered 
at  32°N  and  142°E.  Two  large  areas  of  warm  anomaly  develop¬ 
ment  are  present.  The  first  lies  between  25°-40°N  and  between 
150°-1703E.  The  second  area  lies  off  the  west  coast  of 
North  America  and  extends  westward  along  the  Aleutians  into 
the  Northwest  Pacific. 

3 

Figure  52  shows  the  u*  anomaly  calculated  from  the 

filtered  data  and  integrated  over  the  six-month  period 

(right  side  of  ( I V - 2 )  with  X  =  u  *  )  .  Positive  values  are 

shaded,  indicating  greater  than  normal  storminess  during  the 

period  while  negative  areas  are  unshaded.  The  broad  band  of 
3 

positive  u*  anomaly  values  lies  between  20°-40°N  and  150°E- 
150°W  with  an  extension  of  positive  values  northward  into 
the  Gulf  of  Alaska.  This  pattern  of  storminess  matches  well 
with  that  described  by  Namias  (1978)  for  the  same  period. 
Secondary  maxima  lie  off  the  west  coast  of  Mexico,  in  the 
tropics,  and  in  the  extreme  western  Pacific  (130°E).  Figure 
53  shows  curlzx  anomalies  calculated  from  the  filtered  data 
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and  integrated  over  the  six-month  period  (right  side  of 

( I V  -  2 )  with  X  =  curl^).  The  shaded  areas  indicate  positive 

anomalies  of  curl^  due  to  storminess.  It  is  quite  apparent 

3 

that  the  pattern  is  much  less  coherent  than  the  u*  anomaly 
in  Fig.  52  and  the  general  spatial  scale  is  quite  different 
from  the  SST  anomaly  change  in  Fig.  51. 

In  comparing  Figs.  52  and  53  to  Fig.  51,  the  best 

3 

correlation  can  be  made  between  SST  and  u*  ,  although  there 
is  also  some  relation  with  curl^.  The  large  area  of  SST 
anomaly  decrease  in  the  mid-latitudes  compares  favorably  to 
the  greater  than  normal  amount  of  storminess  in  that  area 

3 

as  measured  by  the  filtered  u*  anomaly.  Also,  the  large 
warm  SST  anomaly  development  in  the  Northwest  Pacific  coin¬ 
cides  with  the  lack  of  storminess  in  that  area  as  indicated 

3 

by  the  negative  u*  anomaly.  A  weak  correlation  also  exists 

3 

between  u*  and  the  cold  SST  anomaly  development  off  the 
California  and  Mexican  coasts,  with  a  better  match  in  the 
West  Pacific  (130°E).  The  area  of  obvious  conflict  lies  in 

3 

the  Gulf  of  Alaska  where  large  positive  u*  values  occur  in 
the  same  place  where  positive  SST  anomalies  developed.  There 
is  a  good  correlation  with  a  negative  anomalous  wind  stress 
curl  in  this  area  but  the  results  are  far  from  conclusive. 
Possible  explanations  for  this  warm  SST  development  may  be 
the  low  frequency  part  of  the  wind  spectrum,  anomalous  surface 
heat  fluxes  or  the  northward  advection  of  warm  surface  water. 

The  purpose  of  this  case  study  was  to  see  if  a 
qualitative  relation  exists  between  SST  anomaly  changes 


and  anomalous  storm  activity.  A  study  is  now  in  progress 
to  determine  the  statistical  relationship  between  monthly 
SST  anomalies  and  synoptic  storm  activity  over  the  entire 
ten-year  period. 


V.  CONCLUSIONS 


The  first  conclusion  drawn  from  climatological  values 

3  3 

of  u*  and  cur1zi  computations  is  that  the  monthly  mean  u* 

calculated  from  filtered  wind  components  is  a  good  indicator 

of  monthly  storminess  and  storm  tracks.  Coherent  patterns 
3 

of  u*  and  tracks  of  associated  synoptic  storms  are  easily 

3 

discernible.  In  addition,  the  u*  computed  from  filtered 

3 

data  constitutes  a  significant  part  (32%)  of  u*.  calculated 
from  the  total  wind  components.  Therefore,  a  major  conclu¬ 
sion  is  that  the  synoptic  storm  part  of  the  frequency  spectrum 
makes  a  large  contribution  to  the  atmospheric  forcing  of  the 
upper  ocean.  On  the  other  hand,  maps  of  climatological  curl^ 
calculated  from  the  filtered  wind  components  are  much  less 
significant  (3.5%),  more  "noisy",  less  coherent,  and  are  not 
as  easy  to  associate  with  storm  tracks.  The  atmospheric  wind 
stress  curl  associated  with  the  more  permanent  events  having 
a  time  scale  greater  than  that  associated  with  synoptic  storms 
appears  to  be  more  important. 

A  second  more  qualitative  conclusion  of  this  study  indi- 

3 

cates  a  relationship  between  anomalous  u*  from  filtered 
data  and  SST  anomalies.  In  the  particular  six-month  period 
of  study,  areas  of  greater  than  normal  synoptic  storm  activity 
were  related  to  the  development  of  negative  or  colder  SST 
anomalies.  Conversely,  areas  of  lower  than  normal  synoptic 
storm  activity  were  related  to  warmer  than  normal  SST  anomaly 


development.  When  comparing  curl2  calculated  from  the  fil¬ 
tered  data  to  the  SST  anomaly  development,  a  much  weaker 
relationship,  if  any,  was  apparent.  It  is  hoped  that  the 
statistical  study  in  progress  will  determine,  more  quantita¬ 
tively,  whether  a  relationship  exists  between  monthly  mean 

3 

anomalies  of  SST  and  u*  calculated  from  the  filtered  data. 

As  noted  earlier,  the  fields  of  curlz  calculated  from 
both  filtered  and  unfiltered  data  have  a  considerable  amount 
of  spatial  noise.  The  noise  comes  from  random  data  errors  in 
time  and  space.  One  solution  to  this  problem  would  be  to 
have  more  accurate  data,  thereby  reducing  the  error  factor. 
Another  solution  would  be  to  increase  the  data  resolution. 

An  increased  resolution  without  increasing  the  accuracy  in 
wind  observations,  however,  would  only  make  the  noise  problem 
associated  with  curl  calculations  worse  or  at  best  unchanged. 
Space  smoothing  could  also  be  used.  But  according  to  Saunders 
(1975),  coarse  data  resolution  (i.e.,  space  smoothed)  may 
underestimate  maxima  and  minima  in  the  curlz  by  as  much  as 
50  percent.  The  ideal  situation  would  be  to  increase  both 
the  accuracy  of  the  data  and  the  space  resolution. 

3 

As  mentioned  earlier,  u*  calculations  from  u  and  v  com¬ 
ponents  included  the  high  pass  effects  only.  This  was  an 
essential  step  in  discussing  storms  and  the  climatology  of 
storms.  Interaction  terms  between  low  and  high  pass  u  and  v 
components  were  not  studied.  The  next  step  is  to  include 
the  interaction  terms  of  low  and  high  pass  filtered  u  and  v 
components  and  then  to  determine  their  effects. 
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Figure  2 
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Response  of  the  41-point  symmetric  high- 
pass  filter  as  a  function  of  frequency. 
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or  February 


Figure  7.  Same  as  Fig.  3  except  contour  values  are 
0 . 5  , 1 . 0 , 1  .  5  , 2 . 0 , 3 . 0 , 4 . 0  .  .  .  .  x  10^  (cm/sec) 
and  the  month  is  May. 


Figure  10.  Same  as  Fig,  7  except  for  August 


MAY  u*  FROM  FILTERED  WINDS 


Figure  19.  Same  as  Fig.  15  except  contour  values  are 

0.5, 1.0, 1.5, 2. 0,2. 5, 3. 0,4. 0,5.0 _ x  10* 

Ccm/sec}3  and  the  month  is  May. 
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Figure  20.  Same  as  Fig.  19  except  for  June. 


NOV  u*3  FROM  FILTERED  WINDS 


Figure  25.  Same  as  Fig.  15  except  for  November. 


DEC  u*3  FROM  FILTERED  WINDS 


Figure  26.  Same  as  Fig.  15  except  for  December. 
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Figure  33.  Same  as  Fig.  31  except  for  July 
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Figure 


35.  Same  as  Fig.  31  except  for  September 
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MAY  CURLZT  FROM  FILTERED  WIND 


Figure  43.  Same  as  Fig.  39  except  contour  values  are 
-4. 0,-3. 5,-3. 0,-2. 5, -2. 0,-1. 5, -1.0, -0.5, 
0.0, 0.5, 1.0, 1.5, 2. 0,2. 5, 3. 0,3. 5, 4.0  x  10" 
dynes/cm3  and  the  month  is  May. 
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Figure  47.  Same  as  Fig.  43  except  for  September 


Figure  48.  Same  as  Fig,  39  except  for  October 
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SST  ANOMALY  CHANGE  FROM  SEP  76  TO  MAR  77 


Figure  51.  SST  anomaly  change  from  September  1976  to 
March  1977.  Contour  intervals  are  0.5°C. 
Shaded  areas  indicate  negative  SST  anomaly 
development. 
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u*  ANOMALIES  FOR  SEP  76  TO  MAR  77 


Figure  52.  Filtered  u*3  anomalies  averaged  from  mid- 
September  1976  to  mid-March  1977.  Contour 
intervals  are  0.5  x  104  (cm/sec)3.  Shaded 
areas  indicate  above  normal  storminess. 


62 


.iiunpiui  nprn 


wpni 


CURL^  ANOMALIES  FOR  SEP  76  TO  MAR  77 


140 


uil  60 


18  0  p  p  pi  60  rfP 

o  o  o  oxd 


1  40  p 


120 


Figure  53.  Filtered  curlzx  anomalies  averaged  over  the 
period  mi d- September  1976  to  mid-March  1977 
Contour  intervals  are  0.5  x  10“9  dynes/cm3. 
Shaded  areas  indicate  above  normal  curlzx 
values. 
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